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Abstract. We present numerical simulations of the hydrodynamical interactions that produce 
circumstellar shells. These simulations include several scenarios, such as wind- wind interaction 
and wind-ISM collisions. In our calculations we have taken into account the presence of dust in 
the stellar wind. Our results show that, while small dust grains tend to be strongly coupled to 
the gas, large dust grains are only weakly coupled. As a result, the distribution of the large dust 
grains is not representative of the gas distribution. Combining these results with observations 
may give us a new way of validating hydrodynamical models of the circumstellar medium. 
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1. Intoduction 

The winds of cool stars contain a large amount of dust. This dust is crucial in driving 
the wind of the star and, at larger distances, it moves along with the flow of the wind 
(Lamers & Cassinelli 1999). When the gas collides with the surrounding medium (inter- 
stellar medium, a previously ejected wind, wind from an other star, etc.) it slows down to 
form a shell, creating a velocity difference between gas and dust. Whether the dust will 
follow the change in velocity depends on the strength of the coupling between gas and 
dust, which in turn depends on such parameters as gas density, velocity difference, tem- 
perature and the size of the dust grains. The larger the dust grains, the more momentum 
they have and the more force will be required to slow them down. In this paper we use 
numerical models of a stellar wind colliding with its surrounding medium to investigate 
the behavior of the dust grains in the resulting shell. 

2. Model 

We use the mpi-amrvac code (Keppens et al. 2011) for our calculations, which solves 
the conservation equations of hydrodynamics on a fixed grid. For these simulations we 
also include the effect of optically thin radiative cooling. The presence of dust is included 
using a two-fluid approximation, where the dust is treated as a gas without internal 
pressure (van Marie et al. 2011), with the interaction between gas described as a drag 
force (Kwok 1975): 



with nd the dust particle density, p the gas density, ad the radius of the dust grains, 
Av the velocity difference between gas and dust and the thermal speed of the gas. 
We simulate the interaction between a stellar wind (10~^ M0/yr, 15km/s) with the 
interstellar medium (ISM), which has a density of 2 particles per cm^. The computational 
domain has a radius of 0.5 pc and an opening angle of 180°. The adaptive mesh refinement 
provides a maximum effective resolution of 3200x800 gridpoints. We assume that the ISM 
contains no dust and that 0.5% of the mass of the wind is in the form of dust grains. 
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Figure 1. Gas mass density in g/cm^ (right panels) and dust particle density (left panels) of 

the circumstellar medium after 10 000 years. The small (0.005 /im) dust grains (left figure) 
show a shell at the contact discontinuity between shocked wind and shocked ISM, indicating 
that the drag force effectively couples them to the gas. The large (0.05 /xm) grains (right figure) 
show no shell, since their larger momentum makes them more difficult to slow down. 

The dust grains are assumed to be made of carbon (internal density 2.9g/cm^). We run 
the simulation twice, once with a grain radius of 0.005 /im and once with a grain radius 
of 0.05 /im. 

3. Results 

Figure [U shows the result of our simulations after 10 000 years. As the gas hits the ISM, 
it creates a shell of shocked gas (shocked ISM on the outside and shocked wind material 
on the inside). The shocked wind, which has the highest density cools very effectively 
(raditive cooling depends on the gas density squared), causing the inner part of the shell 
to be compressed. At this stage the shell does not show instabilities. These will start to 
appear as the shell reaches a larger radius, where the density of the wind becomes less 
than the density of the interstellar medium it sweeps up, making the shell vulnerable to 
Ray leigh- Taylor instabilities. The small (0.005 /im) dust grains show a shell-like structure 
around the contact discontinuity between shocked wind and shocked ISM, which proves 
that the drag force couples them to the gas and has effectively reduced their speed to the 
same velocity as the shocked gas. The larger (0.05 /im) dust grains do not form a shell. 
Their momentum allows them to move almost freely through the shocked gas. 

Whether dust grain distribution will follow the morphology of the circumstellar gas 
depends on the effectiveness of the coupling between dust and gas through the drag 
force. For small grains, this coupling is very effective, forcing the grains to follow the gas. 
For large grains, the coupling is too weak. It is possible to differentiate between grain 
sizes in observations. This will provide us with a new way of validating hydrodynamical 
simulations of the circumstellar environment. 
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